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Inflationary models of the early universe provide a natural mechanism for the 
formation of large scale structure. This success brings to forefront the question of 
naturalness: Does a sufficiently long slow roll inflation occur generically or does it 
require a careful fine tuning of initial parameters? In recent years there has been 
considerable controversy on this issue [1-4]. In particular, for a quadratic potential, 
Kofman, Linde and Mukhanov [2] have argued that the probability of inflation with 
at least 65 e-foldings is close to one, while Gibbons and Turok [4] have argued that 
this probability is suppressed by a factor of ~ 10~^^. We first clarify that such 
dramatically different predictions can arise because the required measure on the 
space of solutions is intrinsically ambiguous in general relativity. We then show that 
this ambiguity can be naturally resolved in loop quantum cosmology (LQC) because 
the big bang is replaced by a big bounce and the bounce surface can be used to 
introduce the structure necessary to specify a satisfactory measure. 

The second goal of the paper is to present a detailed analysis of the infiationary 
dynamics of LQC using analytical and numerical methods. By combining this infor- 
mation with the measure on the space of solutions, we address a sharper question 
than those investigated in [2, 4, 5]: What is the probability of a sufficiently long slow 
roll infiation which is compatible with the seven year WMAP data? We show that 
the probability is very close to 1. 

The material is so organized that cosmologists who may be more interested in the 
inflationary dynamics in LQC than in the subtleties associated with measures can 
skip that material without loss of continuity. 

Key Words: Loop quantum gravity, cosmology, inflation, measures, probability, 
WMAP 7 year data. 



I. INTRODUCTION 

The inflationary paradigm provides a natural mechanism of generating the seeds of inho- 
mogeneities in the cosmic microwave background (CMB) which then evolve to the observed, 
large scale structure of the universe. The general scenario involves a rather small set of 
assumptions: i) Sometime in its early history, the universe underwent a phase of rapid ex- 
pansion during which the Hubble parameter was nearly constant; ii) During this phase, the 
universe was well described by a Friedmann, Lemaitre, Robertson, Walker (FLRW) solution 
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to Einstein's equations together with smaU inhomogeneities which are well approximated by 
first order perturbations; iii) Consider the co-moving Fourier mode ko of perturbations which 
has just re-entered the Hubble radius now. A few e-foldings before the time t{ko) at which 
ko exited the Hubble radius during inflation, Fourier modes of quantum fields describing 
perturbations were in the Bunch-Davis vacuum for co-moving wave numbers in the range 
{ko, ~ 200A;o); and, iv) Soon after a mode exited the Hubble radius, its quantum fluctuation 
can be regarded as a classical perturbation and evolved via linearized Einstein's equations. 
Analysis of these perturbations implies that there must be tiny inhomogeneities at the last 
scattering surface whose detailed features have now been seen in the CMB. Furthermore, 
time evolution of these tiny inhomogeneities produces large scale structures which are in 
excellent qualitative agreement with observations. Therefore, even though the assumptions 
have ad-hoc elements, they appear to capture a germ of truth, not unlike the Bohr atom 
did a hundred years ago. 

For deflniteness, let us assume a quadratic potential for the inflaton and supplement 
our calculations with values of two parameters provided by the seven year WMAP data 
[6]: the amplitude A{t{ki,)) of the scalar power spectrum A^{t{ki,)) and the scalar spectral 
index ns(t(/c*)) at the time the flducial mode A;^ used by WMAP exits the Hubble radius 
{ki, S.5Sko). These numbers, together with the Friedmann equation, determine the values 
of slow roll parameters and the initial data for the inflaton and the gravitational field at time 
t{k^) to within posted observational errors. A slow roll follows and dynamics of perturbations 
during this epoch directly lead to the spectrum of inhomogeneities seen in the CMB. We 
will refer to this inflationary phase as the desired slow roll to distinguish it from other 
inflationary phases that may have occurred, e.g., in an even earlier phase. 

The striking success of the scenario brings to forefront an old issue in a sharper form: Does 
the desired slow roll inflation occur generically in a given theoretical paradigm? That is, do 
generic dynamical trajectories pass through the neighborhood of the values of the inflaton 
and gravitational flelds selected by the WMAP data with its error bars? This would require 
that the inflaton must have been significantly high up compared to the minimum of the 
potential at the onset of the desired slow roll. How did it get there? Is it essential to invoke 
some rare quantum fiuctuations to account for the required initial conditions because the a 
priori probability for their occurrence is low? Or, is the desired slow roll infiation robust in 
the sense that it is realized in 'almost all' dynamical trajectories of the given theory? 

To make these questions precise one needs a well-defined framework to calculate proba- 
bilities of various occurrences within any given theory. A mathematically natural strategy to 
achieve this goal was introduced over two decades ago (see e.g. [7-9]). Recall first that the 
space § of solutions to physically interesting classical systems generally carries the natural 
Liouville measure d/iL. The idea was to calculate a priori probabilities using a fiat probabil- 
ity distribution P{s) = 1 in conjunction with d/i-L- More precisely, the a priori probability 
of an event E is given by the fractional Liouville volume of § occupied by the region R{E) 
consisting of solutions on which the event E is realized [7]. In our case, then, the a pri- 
ori probability is given by the fractional volume occupied by the sub-space of solutions in 
which the desired slow roll infiation occurs. Note that this a priori probability provides 
only a 'bare' estimate and further physical input can and should be used to provide sharper 
probability distributions P{s) and a more reliable likelihood. However, a priori probabilities 
themselves can be directly useful if they are very low or very high. In these cases, it would 
be an especially heavy burden on the fundamental theory to come up with the physical input 
that significantly alters them. 
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However, there is a conceptual obstacle in this calculation: the total Liouville measure of 
the space S of solutions is infinite [9] , hence there is an intrinsic ambiguity in the calculation 
of relative probabilities [3]. But in the observationally favored k=0 FLRW model, this 
divergence is a gauge artifact. More precisely, a gauge group Q acts on S and, although the 
quotient Ei/Q — the space of physically distinct solutions — is compact (with boundary), the 
gauge orbits are non-compact, making the total volume of S infinite. It would first appear 
that the obvious way to avoid the infinite volume is to work directly with the space E>/Q 
of physically distinct solutions. However, as we will see, the Liouville measure does not 
naturally project down to S/Q. Therefore, to calculate probabilities, one has to introduce 
an additional structure. Because the subtleties associated with the interplay between the 
action of the gauge group and the Liouville measure were not well-understood, the necessity 
and importance of this additional step was, apparently, not appreciated. We will see in 
section HIB that there is an intrinsic ambiguity in carrying out this step within general 
relativity. As a recent analysis of Corichi and Karami [10] shows, this ambiguity is directly 
related to the diverging conclusions on probability of inflation in general relativity drawn 
by Koffmann, Linde and Mukhanov [2] and Gibbons and Turok [4]. 

Loop quantum cosmology (LQC) provides a new arena to analyze this issue because 
the big bang singularity is naturally resolved and replaced by a big bounce due to quantum 
geometry effects [11-15]. We will see in section III that, thanks to the presence of a canonical 
bounce time, one can now naturally resolve the ambiguity in the construction of the measure 
on S with finite total volume. Were we to try to mimic this construction in general relativity, 
we would be led to work at the singularity in place of the bounce, where the calculation would 
be meaningless. Away from the Planck regime, LQC is virtually indistinguishable from 
general relativity. However, in the Planck regime, there are huge differences and these are 
crucial in overcoming the obstacle. With this measure at hand, we can calculate the a priori 
probability of the slow roll of the desired type. LQC dynamics are such that this probability 
turns out to be greater than 0.999997: Dynamical trajectories starting from almost all initial 
data at the bounce surface pass through the phase space region selected by the WMAP data. 
Therefore, extreme fine tuning would be necessary to zero-in on solutions where the desired 
slow roll does not occur. 

Some of the results of our investigation were reported in a Letter [5] . Therefore there is 
an inevitable overlap with [5] but there are also key differences. First, in [5] we analyzed 
the likelihood of the occurrence of a slow roll inflation with at least ~ 67 e-foldings in 
the history of the universe to the future of the bounce. In this paper we analyze a much 
sharper question: What is the probability of occurrence of a slow roll with initial conditions 
that are compatible with the 7 year WMAP data? Thus we now focus only on that slow 
roll phase which is directly relevant to structure formation. Second, numerical simulations 
reported in [5] used values of cosmological parameters — the mass of the inflaton and the 
values of the slow roll parameters — from Linde's 2006 review [16] while in this paper we 
use instead the more recent results of the 7 year WMAP data [6] . This accounts for some 
differences in some of the detailed numerical results. Finally, and more importantly, our goal 
now is broader than that of [5] in the following sense. In LQC, the big bounce is followed 
by a qualitatively new phase of super-inflation which could have observable consequences. 
Therefore, it is important to have a sufficiently detailed account of the new dynamics from 
the big bounce to the onset of the desired slow roll. A second goal of this paper is to provide 
this analysis. This detailed description is likely to serve as the point of departure of further 
work bridging the Planck era of LQC to the inflationary paradigm described in the beginning 
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of this section. This bridge may, for example, provide a better understanding of why the 
quantum state of relevant modes is well approximated by the Bunch Davis vacuum at the 
onset of inflation, and may even furnish the quantum gravity corrections to this state [17]. 

The paper is organized as follows. In section II we first recall the relevant features of LQC 
and then introduce the phase space and basic equations. In section III we first introduce 
the Liouville measure on the space S of solutions, discuss the issue of gauge and obtain the 
measure on the space S/Q of physically distinct solutions in LQC. In section IV we discuss 
in detail the LQC dynamics from the big bounce to the end of infiation using a combination 
of analytical and numerical methods. Using this information and the measure introduced in 
section III, we calculate the a priori probability of obtaining the desired slow roll. As noted 
already, in LQC this probability is very close to 1. Section V summarizes the main results 
and compares and contrasts them with related results in the literature. 

The material is organized so that cosmologists who may be more interested in infiationary 
dynamics of LQC than in the issue of measures can skip sections II B and especially III 
without loss of continuity. 

II. PRELIMINARIES 

This section is divided into two parts: In the first, we recall the distinguishing features 
of LQC that are important to our analysis. In the second, we present effective LQC equa- 
tions governing dynamics of the FLRW model coupled to a scalar field with any potential 
(satisfying mild regularity conditions) and introduce the Liouville measure on the space of 
solutions. 

A. Distinguishing features of LQC 

In LQC, one apphes the basic principles of loop quantum gravity (LQC) [18-20] to simple 
cosmological models. Thanks to the quantum geometry underlying LQC, LQC differs from 
the older Wheeler-DeWitt theory already at the kinematical level. It turns out that the 
Wheeler-DeWitt equation is no longer well-defined on the new kinematical Hilbert space. 
Instead, now the quantum Hamiltonian constraint has to be obtained via a procedure that 
pays due attention to the quantum geometry of LQC, in particular, the area gap A = 
TT^fpj. (Here 7 is the Barbero-Immirzi parameter of LQC, whose value 7 ~ 0.24 is 
fixed by black hole entropy calculations [21, 22].) Somewhat surprisingly, the resulting 
dynamics naturally resolves the big bang and big crunch singularities of general relativity 
[23]. Exotic matter is not needed; indeed matter fields can satisfy all the standard energy 
conditions. Detailed analysis has been carried out in a variety of models: the k=0, 1 FLRW 
space-times with or without a cosmological constant [12-14]; Bianchi models [24-26] which 
admit anisotropics as well as gravitational waves; and Gowdy models [27] which admit 
inhomogeneities, and therefore an infinite number of degrees of freedom. The FLRW models 
have been studied most extensively, using both analytical and numerical methods to solve 
the exact quantum equations of LQC [12, 13, 15]. In these models, the big bang and the 
big-crunch are replaced by a quantum bounce, which is followed by a robust phase of super- 
infiation. Interestingly, full quantum dynamics, including the bounce, is well-approximated 
by certain effective equations [12, 15, 28]. These equations imply that all strong curvature 
singularities — including the big rip and sudden-death — are resolved in FLRW models with 
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matter satisfying an equation of state of the type p = wp [29]. (For recent reviews, sec [30].) 

In this paper we will restrict our matter source to be a scalar field with the standard 
(positive) kinetic energy and a suitable potential. Since all the prior discussion of probabil- 
ities is based on general relativity, to facilitate comparison we use effective equations rather 
than the full quantum theory. Finally, we will use the natural Planck units c=h=G=l 
(rather than 87rG=l, often employed in cosmology). Planck length will be denoted by ip\ 
and Planck mass by mpi. The fundamental time unit, Spi := ^/ Gh/c^, will be referred to as 
a Planck second. 

In LQC, it is convenient to encode spatial geometry in a variable v proportional to 
the physical volume of a fixed, fiducial, cubical cell, rather than the scale factor a. The 
conjugate momentum is denoted by b. These are related to the scale factor and its conjugate 
momentum via^ 



and 6=-^--!^ (2.1) 



and 



where Vq is the co- moving volume of the fiducial cell, so that its physical volume is o^Vq. 
(Thus, the only non-vanishing Poisson bracket is {v, b} = —2.) On solutions to Einstein's 
equations, b is related to the standard Hubble parameter H = a/a via b = [15]. However, 
LQC modifies Einstein dynamics and on solutions to the LQC effective equations we have 

H = ^ sin2A& ^ (0.93 mpi) sin2A6 (2.2) 

where := A ^ 5.2£pj is the 'area- gap' that sets the discreteness scale of LQC. b ranges 
over (0,7r/A) in LQC and general relativity is recovered in the limit A — ?■ 0. 

Quantum geometry effects modify the geometric, left side of Einstein's equations. In 
particular, the Priedmann equation becomes 

sin^ A6 Stt 8tt ,6'^ ,x\ „x 

To compare with the standard Priedmann equation — (Stt/S) p, it is often convenient to 
use (2.2) to write (2.3) as 

where pcrit = 3/87r7^A^ ~ 0.41ppi. By inspection it is clear from Eqs (2.2) - (2.4) that, away 
from the Planck regime — i.e., when A6 <S 1, or, p <C Pcrit — we recover classical general 
relativity. However, modifications in the Planck regime are drastic. The main features of 
this new physics can be summarized as follows. 

• In general relativity, the Priedmann equation implies that if the matter density is pos- 
itive, a cannot vanish. Therefore every solution represents either a contracting universe or 



In LQG, the basic variable is a triad rather than a 3-metric and in the LQC hteraturc v is taken to be the 
oriented volume which is positive for positively oriented triads and negative for negatively oriented ones. 
However, since the change of orientation is a large gauge transformation, in the classical and effective 
theories on can restrict oneself just to positive v. We have done so for simplicity of discussion. 
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an expanding one. By contrast, the LQC modified Fricdmann equation (2.4) impfies that 
i) vanishes at p = Pcrit- This is the quantum bounce. To its past, the solution represents a 
contracting universe with v < and to its future, an expanding one with v > 0. 

• As is customary in the hterature on probabihties, let us ignore the exceptional de 
Sitter solutions. On all other solutions h is monotonically non-increasing, evolving from 
from 6 = tt/A in the infinite past to in the infinite future. Eqs (2.3) and (2.4) imply that 
6 = 7r/2A at the bounce. Thus, each solution undergoes precisely one bounce. 

• In contrast to general relativity, the Hubble parameter H — v/2)V\sno longer monotonic 
in LQC. It vanishes at the bounce while in general relativity it diverges at the singularity 
and is large in the entire Planck regime. In LQC, H is bounded above, \H\ < 0.93 mpi, and 
achieves its upper bound in every solution at the end of super- inflation. 

• If the potential V{(p) is bounded below, say V > V^, then it follows from (2.3) that 0^ 
is bounded by 2pcrit — '^Vo- H V grows unboundedly for large then |0| is also bounded. 
For example, for V — m^0^/2, we have m|0|max = 0.90 m|j. 

• When the potential is bounded below, \H\ is bounded above by 10.29 mpj. The Ricci 
scalar — the only non-trivial curvature scalar in these models — is bounded above by 31mp[. 
Thus, physical quantities which diverge at the big bang of general relativity cannot exceed 
certain finite, maximum values in LQC. One can also show that ii v ^ initially, it cannot 
vanish in finite proper time along any solution. Thus, the LQC solutions are everywhere 
regular irrespective of whether one focuses on matter density, curvature or the scale factor. 



B. The Liouville measure 

The full set of space-time equations of motion can be written in terms of v{t), (f){t). These 
variables are subject to the constraint (2.4) and evolve via: 

v^^[{p- V{<t^)f + V(0)(pcrit - V{<p))\ (2.5) 

Pcrit 

0+^0 + y^ = O. (2.6) 

To calculate probabilities, in section III we have to equip the space S of solutions to these 
equations with a natural measure. As a first step in that procedure, we will now obtain 
a phase space formulation of these equations. The phase space and equations of general 
relativity can be recovered by taking the limit A — )■ (which in particular implies pcrit ~^ oo). 

The phase space V consists of quadruplets (f , 6; (f),P(j)), where is given by; = 
2TTjipiV 4>. The variables (p,Pij) range over the entire real line, v over the positive half 
of the real line, while b e [0,7r/2A] (since we focus only on the post-bounce branches of 
solutions). Thus, the symplectic 2- form is given by 

fl = d(j)Adp^ + -db A dv (2.7) 

Hence the Liouville measure on F is simply 2d/iL = d^dp.^dfcdi'. 

The LQC Friedmann equation implies that these variables must lie on a constraint surface 
F defined by 

C= --4^ sin^Afc + -^ + 2777^^(0) 0. (2.8) 

47 A2 ATT-fV 
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They evolve via 

0=^, Pi^^ = -27rjvV^. (2.10) 

As is well-known, the space of solutions § is naturally isomorphic to a gauge-fixed surface, 
i.e., a 2-dimensional surface F of F which is intersected by each dynamical trajectory once 
and only once. Since b is monotonic in each solution, an obvious strategy is to choose for F 
a 2-dimensional surface b — bp {b: fixed constant) within F. 

It is straightforward to pull-back the symplcctic structure to this 2-dimensional gauge- 
fixed surface F. Using the constraint (2.8), it is convenient to coordinatize F using v, (j) and 
express the pulled-back symplectic structure in terms of them: 

Q77- 1 

Q = [— sin^ Xbo - 8TT^-fW{(f))] ^ d(f)Adv (2.11) 
A 

This 2-form provides provides a Liouville measure d/(L on F and hence of the space § of 
solutions to the effective equations, given simply by 

d/iL = sin^ Xbo - 87r^-f^V{<j))] ^ d<f)dv. (2.12) 

The most natural choice in LQC is to set bg = 7r/2A so that F is just the 'bounce surface'. We 
will make this choice because it also turns out to be convenient for calculations. However, 
since the dynamical fiow preserves the symplectic structure on F, this measure on § is 
insensitive to the choice of bo used in gauge fixing. 



III. PROBABILITY CONSIDERATIONS 

Recall from section I that the a priori probability of occurrence of any event E is to be 
given by the fractional volume of the region R{E) in § spanned by solutions in which E 
occurs: 

P^E^^lm)^ (3.1) 

Since the Liouville measure is purely kinematical, this method of calculating a priori prob- 
abilities realizes Laplace's principle of indifference [31]. This interpretation can be made 
explicit as follows. Physical input can provide a (non-negative) probability density function 
p(s) on S satisfying the normalization condition [J^pdfii^]/[J^dfii] = 1. The corresponding 
(more refined) probability of occurrence of E is then given by 

fF/(F\ Pi-^) ^Al 

Pp{E) = '7y^^. . (3.2) 

Now, one can quantify the information contained in p{s) via 

_ /gp(g) Inp(g) dfiL 

/,p(.)dAL ■ ^'-'^ 
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Since / is minimized by p(s) = 1, (3.1) is the a priori probability, free of any input or bias. 
As mentioned in section I, these 'bare' probabihties are useful when they are extremely low 
or extremely high. In these cases it is a heavy burden on any theory to provide sufficient 
information to significantly overcome this bare probability. 

In section III A we illustrate these ideas using a 2-dimensional simple harmonic oscillator 
constrained to have a fixed energy. In this case the solution space S is compact. Since the 
total Liouville volume of § is finite, one can directly implement the ideas outlined above 
to calculate various probabilities. In section IIIB we turn to LQC. In this case, S is non- 
compact. But we will see that this non-compactness can be directly attributed to the action 
of a gauge group Q on §. Therefore, it is now natural to work with the space S S/Q of 
physically distinct solutions which is compact. However, there is a subtlety which makes the 
calculation of desired probabilities ambiguous in general relativity. We discuss this problem 
and show that it has a natural resolution in LQC. 

Cosmologists who are more interested in the LQC dynamics than in the issue of measure 
on § can skip section III without loss of continuity. 



A. Harmonic oscillator 

Since we encounter a constrained Hamiltonian system both in general relativity and (ef- 
fective) LQC, let us begin with a simple example with this feature to illustrate the procedure 
of calculating probabilities. Consider a 2-dimensional simple harmonic oscillator constrained 
to have a fixed energy E. Let us fix the mass and spring constant to unity for simplicity. 
The topology of our phase- space F is then W^, in which we can choose canonical coordinates 
Xi,Pu X2,P2, with symplectic structure: 

uj — dxi A dpi -|- dx2 A dp2 (3.4) 

The Hamiltonian constraint is: 

H:=^{xI+pI + xI+pI) = E (3.5) 

The constrained surface T is the 3-sphcre with radius \/2E in and thus compact. The 
Hamiltonian vector field is given by 

d d d d _ „. 

-^H^Xi- pi- \-X2^ P2^— (3-6) 

Opi OXi Op2 0X2 

The vector field is of course tangential to V and its orbits provide a Hopf fibration: each 
orbit is an fiber with the base space being S"^. This base space represents the space S of 
solutions of the constrained system under consideration. 

It is convenient to use a set of coordinates adapted to this Hopf fibration: 

Xi -\- ipi = \PlE e*^^ sin r\ = z\ and X2 + ip2 = y/2E e*^^ cos rj = Z2 (3.7) 

where rj G (0,7r/2) and ^1,^2 G (0,27r). Then T is given by l^ip -|- \z2\'^ = 2E and the angles 
^7)^1)^2 provide a a natural set of intrinsic coordinates on it. The pull-back of to f can 
now be expressed as 

n^2E sin 2r] dr] A d^" (3.8) 
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and the restriction of the Hamiltonian vector field to F is given by 

X, = 2 A (3.9) 

where we have set 2^=*= = ± ^2- As expected, Q is degenerate and the degenerate direction 
is precisely Xh- It is obvious that ry,^~ and VL are Lie dragged by Xh- Hence they have 
unambiguous projections to the 2-sphere of orbits of Xh, i.e., to the space § of solutions. 
The total Liouville volume of S is given by: 

Cl ^ 2E dr]sm2r] dC = '2n{2E) . (3.10) 
Js Jo Jo 

We can now use the Liouville measure d/iL = (2ii^) sin2r7d77d^^ on S to calculate the 
a priori probabilities of physical events of interest. Consider example, the total angular 
momentum 

J := X1P2 — X2P1 = —E sin2?7 sin2^~ (3-11) 

which is a Dirac observable and thus projects down from F to S unambiguously. Another 
Dirac observable is the eccentricity of the orbit in xi,X2 plane defined in each solution, 
e = Tmin/fmax whcrc t"^ — x\ -\- x\. One can show that e is determined completely by the 
ratio \J\lE: 



1 J2 



(3.12) 



To gain insight into the likelihood of various values for both observables, it is natural to 
ask, e.g., for the a priori probability that \J\/ E is larger than a given number / e [0, 1]. We 
obtain : 

P{(\J\lE)>f)^^jjl^l-f (3.13) 

where A is the region such that sin2?7sin2^~ > /. We find that the probability for e > 1/2 
is P|j|/E>i6/25 = 1 - 16/25 = 9/25 and the probability for e > 3/4 is only 49/625 ^ 0.08. We 
can therefore conclude that almost circular orbits (e ~ 1) are rare according to our measure, 
and thus if the orbit of a particle in this system were observed to be very close to circular, 
further physical explanation would be required. 

Finally, in this analysis we identified the solution space § with the quotient of f by the 
orbits of the Hamiltonian vector field. In general relativity and LQC it is more convenient 
to gauge fix. In the present case, because the orbits of Xh provide a Hopf fibration of F a 
global gauge fixing is not available. Still, for calculating probabilities, we can ignore sets of 
measure zero and parameterize S by points of a 2-dimensional surface F given by ^+ = const. 
For physical questions such as the ones discussed above, the resulting a priori probabilities 
are insensitive to the choice of the constant and the results are the same as those obtained 
above. 



B. Measures for general relativity and LQC 



Recall from section II that the phase space F of LQC is naturally coordinatized by 
{v, b; 0,p^) where v takes values in the positive half line, 0,p<^ on R and 6 e (0, tt/A). These 
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phase space variables are subject to the Hamiltonian constraint (2.8) of effective LQC and 
dynamics of this theory is generated by the Hamiltonian vector field Xq of the constraint 
function C: 



.3v sin2A6. d / . d , . d , t/ ^ ^ ('i^A\ 

There is a gauge symmetry which will play an important role in our probability considera- 
tions. Since v is, by definition, the physical volume of a fiducial cell and the choice of this 
cell is arbitrary, physics does not change under the transformation 

{v,b;(f),p^) — > {av,b;(f),ap^) (3.15) 

which corresponds merely to a rescaling of the fiducial cell by a positive number a (recall: 
P0 = 27r7£pj V 0) . Thus the phase space F is acted upon by the action of a gauge group Q 
generated by (3.15). (This gauge freedom is equivalent to that in rescaling the scale factor 
in the more familiar Wheeler-DeWitt phase space.) At an infinitesimal level this rescaling 
defines a vector field 

Under this rescaling the constraint function (2.8) and the symplectic structure (2.7) are 

merely rescaled, CqC = C and Cg^ = such that the Hamiltonian vector field (3.14) is 
left invariant, CqXc = 0. Hence, as one would expect of gauge transformations, Q descends 
to the space S of solutions. These considerations hold both for effective LQC and general 
relativity. 

Next, let us find a suitable parametrization of §. Note first that we can solve (2.8) for 

P^-^[^ - SttV^ V(0)] (3.17) 

Without loss of generahty we can restrict ourselves to the positive branch > because 
identical considerations will apply to the negative branch. Thus the constraint surface F is 
coordinatized by (v,6, 0). Since the dynamical trajectories lie in f, the space of solutions 
§ is 2-dimensional and naturally isomorphic to any gauge-fixed surface F in F which is 
intersected by each dynamical trajectory once and only once. Since b is monotonic along 
dynamical trajectories, we can take F to be the 2-dimensional surface {v,b — bo,(f)) for any 
fixed constant bg in the domain (0, tt/A) of The Liouville measure d/tL on F, and hence 
on S, is then given by (2.12): 

37r 1 
d/fL = sin^ Xbo - 8Ti^^^V{(j))] ^ d(f>dv. (3.18) 
A 

Since it is induced by the symplectic structure, the d/tL volume of any region on the b — bo 
section is the same as that of its image on any other b = const surface under the Hamiltonian 
flow Xc- In this sense, d/iL is insensitive to the choice of bg- The analogous structures 



^ It is then convenient to choose the freedom in rescahng the constraint i.e., in the choice of the lapse 
function — to make b the afHne parameter of the Hamiltonian vector field, so that Xq now maps each 
b = const surface to another b = const surface. 
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Tgr, f gr, Tgr and d/iL iii general relativity are obtained by taking the limit A — >■ 0. In 
particular, in this limit b H, the Hubble parameter. Hence in general relativity one can 
naturally fix the gauge by setting H = Ho for any positive constant Ho- 

The essential problem in both cases is that the spaces of solutions § and Sgr are non- 
compact and their volume with respect to d/(L and dp,^^ is infinite. Therefore we cannot 
directly use the procedure of section HI A to calculate probabilities of events of interest. Let 
us restrict ourselves to positive potentials V{(f)) — such as the quadratic and quartic ones 
which are often analyzed in detail — which diverge at = ±oo. Then the constraint (2.8) 
(and its GR analog, the Priedmann equation) immediately implies that the range of is 
bounded on F. Thus, the non-compactness of F and Fgr arises only because v ranges over 
the entire positive real line. However, as we saw in the beginning of this sub-section, the 
transformation {v,(j)) {av, 4>) is a gauge motion. Thus, the total Liouville measure of § is 
infinite simply because each orbit of the gauge group Q is infinitely long. 

Let us spell out the interaction between the gauge group and the space of solutions S 
explicitly. Since CqC = 0, the gauge vector field G is tangential to F. Furthermore since 
Cob = (and CqH = in general relativity), G is also tangential to the gauge-fixed surface 
b = bo in effective LQC (and to the gauge- fixed surface H = Ho 'va. general relativity). Thus, 
G is tangential to F (and Fgr), i-e., the action of the gauge group Q naturally respects our 
parametrization of the space § of solutions using F. The space of physically distinct solutions 
is the quotient § := S/^, coordinatized just by which takes values on a compact interval, 
say (0min, 0max)- Siucc wc are interested in gauge invariant questions, regions R{E) of S 
consisting of solutions in which physical events — such as the desired slow roll infiation — 
occur contain whole orbits of Q and project down unambiguously to §. To calculate the 
prol)al)ility of such events, then, it would suffice to have available a natural measure on S. 
Since § is a closed interval, one would expect a natural construction to provide it with a 
finite total measure, making the calculation of desired probabilities well-defined as in section 

niA. 

Thus, the natural strategy is to attempt to introduce a volume element on S = V/Q 
starting from the symplectic structure Ct on F."^ The obvious candidate for such a volume 
element is the natural 1-form tj^ = which is transverse to the gauge direction G". 

Unfortunately, since CcCoa = {i^G^ai3)G^ = cJq 7^ 0, this 1-form Coa docs not project down 
to S. Therefore to implement this strategy one has to introduce an additional structure. Let 
us fix a line So in S = F given hy v — Vo for some constant Voi which is naturally isomorphic 
to^ — T/Q.^ One can then just pull-back the 1-form u to So and perform integrations there. 
Expression (2.11) of Q, imphes that the total volume of So is given by 

r /■(/'max Q„ 1 

/ a;= / (— sin2A6o-87rV^(0))'^^o#, (3-19) 

J So •/— 0max 

which is manifestly finite. Finally, a change in the value of Vo would simply rcscalc Cj by 



^ A straightforward push-forward of the Liouville measure /Ul on § = f is not useful because it would assign 
to any interval onT/Q the same volume as its inverse image in T which is infinite. 

^ The use of a more general line, v = f((j)), would correspond to giving the gauge orbits a (j> dependent 
regularized volume. This would introduce an ad- hoc new input in the construction. With v — Vo, the 
regularized volume of each orbit is the same and, as explained below, it drops out in the calculation of 
probabilities of physical events. 
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a constant. Since the relative probability P{E) of the occurrence of a physical event E is 
given by 

P(E) = , (3.20) 

where is the intersection of the v = Vo line and the region R{E) on which the property 

E is realized, clearly, P{E) is independent of the choice of Vg made above. 

Thus, it would appear that we have satisfactorily eliminated the spurious infinity in the 
total Liouville measure and arrived at a natural procedure to calculate desired probabilities. 
However, there is an important subtlety which makes the procedure ambiguous. Recall 
that in order to arrive at cu we made two gauge choices: We set b = bg (or, H = Hg in 
GR) and then v — Vg. The symplectic structure and the 1-form a) on S is preserved by 
the Hamiltonian flow and (with the appropriately chosen lapse) the flow maps each of our 
2-dimensional gauge fixed surfaces b = const to another such surface. However, under this 
flow, the data {v = Vo,b = bo] (l),p^ = 'P4>{vo, bo, 0)) with constant v on the b = bo slice is not 
mapped to data with constant v at another slice, say b = bi. (Here p^ivo, bo, (p) is the solution 
(3.17) of the constraint). Therefore the final step, where we carry out the integration on 
the V — Vo line, of our procedure sensitively depends on our initial choice of bo (or H = Ho 
in general relativity). The question then is: Is there a canonical choice bo one can make in 
LQC (or. Ho in general relativity)? 

In general relativity, the answer is in the negative: there is no preferred value of the Hub- 
ble parameter or a canonical instant of time in any solution. Therefore the calculation of 
probabilities has an intrinsic ambiguity. Although there are also other differences (discussed 
in section V), this ambiguity lies at the source of dramatically different predictions on the 
probability of inflation because, in effect, by choosing very different values of Ho, one can 
arrive at very different measures for computing probabilities. In terms of these considera- 
tions, Gibbons and Turok chose a low value of Hg and found that the sub-space of solutions 
admitting a sufficiently long slow roll occupies an extremely tiny relative volume of S with 
respect to the resulting measure on S. A higher value of Ho — as was advocated in the 
early literature (see, e.g. [32]) — increases this probability very substantially. In LQG, on 
the other hand, the bounce surface provides a canonical 'time' instant, i.e., a natural value 
for bo, namely bo = 7r/2A where the matter density attains its maximum value. The corre- 
sponding surface in general relativity would be the big bang singularity and, unfortunately, 
we cannot use the 'data at the big bang' to parameterize the space S of solutions. ^ 

To summarize, as discussed in the early literature [7-9], the Liouville measure on the 
phase space naturally descends to the space § of solutions but the total d/tL measure of 
§ is infinite. In the observationally favored k=0 models this infinity is physically spurious 
because it arises only because the length of gauge orbits is infinite. For physical questions 
one can work with the quotient S/Q which is just a bounded interval of the real line both in 
general relativity and LQG. However, although the Hamiltonian vector field Xq is invariant 
under the action of the gauge group Q, the symplectic structure is not. Gonsequently, the 
Liouville measure fails to naturally descend to a measure (with finite total volume) on E>/Q. 
To introduce such a measure, an additional structure is needed. Because of the presence of 



^ A common suggestion that one should carry out the calculation when the Hubble parameter is of Planck 
scale, rather than infinite, comes close to the natural strategy in LQC but the suggestion is only qualitative 
and, strictly speaking, quite different because in LQC the Hubble parameter vanishes at the bounce surface. 
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the bounce, this structure is naturally available in LQC In general relativity, by contrast, 
there does not appear to exist a natural way to introduce the required structure; the analog 
of the bounce surface would be the big bang singularity itself! Consequently, in general 
relativity, the calculation of a priori probabilities of various physical events along the lines 
of [7-9] is intrinsically ambiguous. 

Remark: In LQC, the final result can be stated directly in terms of the Liouville measure 
d/iL on S = r. The regions of interest R{E) spanned by solutions in which a physical event 
E occurs contain whole orbits of the gauge group Q: R{E) — I x where / is a closed 
interval in [— 0max, 0max] and denotes the v axis (0 and v evaluated at the bounce surface 
b = 7r/2A). The associated probability P{E) is given by first introducing a 'slab' /^^ in the 
2- dimensional space F, bounded hy v — Vo and v — 1/vo with < Vo < 1, and then taking 
the limit Vg ^ 0: 

= Liouville Volume of [I x U 

vo-^o Liouville Volume of [/total x ^oJ 
This is the expression that was used in [5]. 



IV. INFLATIONARY DYNAMICS IN LQC 

The event E of interest to us is the presumed slow roll that is compatible with the WMAP 
data. With a canonical measure on the space S/Q of physically distinct solutions of LQC 
at hand, we can now ask for the relative volume occupied by the region R{E) consisting 
of solutions in which E occurs. In this section we will first characterize this subspace by 
analyzing in detail the post-bounce dynamics of LQC and then calculate its relative volume. 

In section IV A we collect useful facts about the desired slow roll that arc used throughput 
the subsequent subsections. Solutions of effective LQC equations can be naturally identified 
with their initial data at the bounce surface. It turns out that the qualitative features 
of the relevant dynamical evolution are largely dictated by the ratio of the kinetic to the 
potential energy in this initial data. In particular, as we will see in section IV C, only those 
dynamical trajectories for which the kinetic energy overwhelmingly dominates fail to lie in 
R{E). To bring out this point, in section IV B we will study in detail the solutions in which 
the kinetic energy completely overwhelms the potential energy at the bounce. In section 
IV C we first briefly discuss dynamics in other solutions and then show that except in the 
case of very extreme kinetic energy domination at the bounce, the solutions necessarily lie 
in R[E). In section IV D we find a bound on the relative volume of R{E) which shows that 
the probabihty of seeing the desired phase of slow roll inflation is extremely close to 1 in 
LQC. While the detailed discussion uses a quadratic potential, we argue that the qualitative 
result is likely to hold much more generally. 



A. The desired slow roll 



Let us begin by recalling the slow roll parameters that are used in the discussion of the 
WMAP data. There are two conceptually distinct sets of parameters [33] 
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Smallness of the first set ensures that the Hubble parameter is changing very slowly in the 
dynamical phase under consideration. (Sometimes one uses the symbol en for e and rju for 
T] to highlight this.) Inflation — i.e., accelerated expansion — is characterized by e < 1. For 
slow roll, on the other hand, we need e and <S 1. As is common in the literature, in 
our study of slow roll, we will keep terms of order 1 in these parameters but ignore higher 
order terms. For concretencss, we will end the slow roll phase when e = 0.1 so that the 
errors will be less than a percent. In theoretical investigations one often uses the second set 
of slow roll parameters which is tailored to the properties of the potential. During the epoch 
in which general relativity is an excellent approximation, the two sets are closely related by 
dynamics. In particular, ey — [(1 + r)/3)/{l — e/3)]^ e. Hence the difference between e and 
ey is of second order in e, rj. Finally, these parameters are best suited for studying dynamics 
during the general relativity era. In full LQC, on the other hand, H = ai the end of super- 
inflation but the super-inflation phase is far from being quiescent if the bounce is kinetic 
energy dominated. A better definition of the first slow roll parameter is e = (3KE)/(2p). 
In general relativity, this definition agrees with (4.1) but with this definition the slow roll 
condition e <C 1 is violated even at the end of super-inflation if the bounce is kinetic energy 
dominated. In this paper we will use the standard definition (4.1) because all our slow roll 
considerations will refer only to the general relativity epoch. 

Next, let us recall the constraints imposed by the WMAP data on scalar perturbations. 
The data is tailored to the co-moving wave number A;* given by [6] 

k 

— = 2 X 10"''^ Mpc"^ or = 8.58 ko (4.2) 

where refers to the scale factor today and, as before, kg refers to the wave number that has 
just re-entered the Hubble radius today. (It is only the combination fc^/oo that has direct 
physical meaning; 2110,0/ ko is the physical wave length of this reference mode today.) Within 
infiationary models, the data constrains initial values of fields describing the homogeneous 
isotropic background at time t(A;*) i.e., the time at which the mode k^, exits the Hubble radius 
during inflation. In the rest of this subsection, we will make these constraints explicit. 

The amplitude A{t{ki,)) of the scalar power spectrum A^(A;^) at this wave number is 
given by: 

and the scalar spectral index ns{ki,) is given by 

ns{t{k,)) = 1 - = 0.968 (4.4) 

dlnA; fc* 

with error bars of about ±4.50% for A and ±1.25% for us [6]. 

In the next two sub-sections, we will focus on the quadratic potential, V{(j)) = (l/2)m^0^. 
The equation of motion of the scalar field (2.5) then simpfifies to: 

+ 3/^0 + = 0. (4.5) 

For this potential, 1 — = 4e (irrespective of the value of the infiaton mass) so we conclude 
that the value of the slow roll parameter at the time t{k^) is given by: 

e(i(A;0) = 8 X 10-^ (4.6) 
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The observed value of the amplitude of the scalar power spectrum then determines the 
Hubble parameter H{t{k^)) and the Hubble radius i?H(^(^*)): 

H(t(h)) = 7.83 X 10-^mpi or Rnitih)) = 1-28 x 10^£pi (4.7) 

Therefore, it follows from our discussion in section HA that the LQC corrections to general 
relativity are highly suppressed at t — and thereafter. For simplicity, we will use this 

fact in what follows. 

Using the approximation e = ey at time t{k^,) (which is consistent with observational error 
bars) we can determine (f){t{ki^)) from (4.6). Then, using the implication e/3 = 0^/(0^ + 
m'^cjp') of the definition of e, and the Friedmann equation, it is easy to obtain the following 
values 

(t){t{k^)) = ±3.15 mpi, ^{t{k^)) = ±1.98 x 10"^m|i 
m = 1.21 X IQ-^mpi ri{t{k^)) = 1.61 x IQ'^ . (4.8) 

(In the reduced Planck units, often used in the cosmological literature, the inflaton mass 
is m = 6.06 X 10~^Mpi.) Because this value of m is somewhat different from that given 
in [16], which was used in [5], some of the details of numerical results in this section differ 
from those reported there. However, as noted in [5], the main results reported there do not 
change appreciably even if this value is changed by a couple of orders of magnitude. 

For the desired slow roll, without loss of generality we can assume (f){t{k^)) < and 
(0)(t(/ov)) > 0, so that decreases as the inflaton slides down the potential. Then we can 
ask for the number of e-foldings of slow roll inflation starting from this time, t = t{k^) when 
e = 8 X 10~^ until it increases to 0.1. This is the 'desired' slow roll phase of interest to our 
analysis. We have: 



0cnd 8^ /--Acnd (1 + 2) y 



M := In = / ^ d0 = — ^ / ' li ^ d0 

a{t{k^)) J 0(t(fcO) ^Pi J mk*)) (1 - 3) ^ 

" ^ + l^^^ + |) - ^ m^k.)) - Cd] - 57.5 (4.9) 

where in the second line we have ignored second and higher order terms in the slow roll 
parameters.^ For simplicity, here we have used the general relativity field equations. The 
LQC result 

gives corrections of order (0/0max)^ ~ 10"-'^° which are totally negligible. 



^ Consistency with assumptions made in the inflationary scenario require the actual slow roll to begin 
somewhat earlier because ki, = 8.58ko and the use of Bunch-Davis vacuum requires the mode kg to be 
well within the Hubble radius at the onset. The initial data we begin with ensures that this will be the 
case. If we assume that at the onset the Rn is 100 times the physical wavelength of the mode kg, then 
there are ~ 6.75 e-foldings between the onset and t{ki,), bringing the total number of slow roll e-foldings 
from the onset until the end (i.e. until e = 0.1) to A/totai ^ 64. This is a lower bound; an earlier onset 
and hence a larger number of e-foldings is permissible within this scenario. 
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To summarize, assuming the inflationary scenario with a quadratic potential, the WMAP 
data provides us within observational errors of < 2% : i) the value of the inflaton mass; 
and ii) the initial data for inflation at some time t{k^) in the early history of the universe. 
This data then automatically leads to the desired slow roll inflation with ~ 57 e-foldings 
starting t — t{ki,). The question now is: What is the probability that the LQC dynamical 
trajectories pass through this very small region of the phase space, irrespective of what their 
initial data are at the bounce surface? 



B. Extreme kinetic energy domination at the bounce 

In this sub-section we will study in detail the LQC dynamics of solutions for which less 
than 10~^° of energy density is in the potential at the bounce. Recall that, since the total 
matter density at the bounce always equals Pcrit, the value |0b| of the scalar field at the 
bounce point is bounded above by 0max = 0.90 rripjm ~ 7.47 x lO^mpi. We will use the 
fraction 




(4.11) 



to bin the solutions. By definition, / takes values in the interval [—1, 1]. Note also that 
= V{(p)/pcnt, the fraction of the total energy density at the bounce that is in the po- 
tential. Therefore in this section we will focus on solutions with |/| < 10~^. Starting at 
the bounce, we will first describe the main features of the LQC evolution using suitable 
analytical approximations and compare these with results obtained by numerical evolution 
of the exact system. 

Let us begin by noting a symmetry of the phase space equations of motion: Given a 
solution (0(t),p</,(t); v{t), b{t)) to (2.8) and (2.9), (— ^(t), —P(j,{t); v{t), b{t)) is also a solution. 
Therefore, in the discussion of dynamics it suffices to focus on the initial data at the bounce 
point where p^\b > 0, i.e., 0b > 0, allowing (f)B to take both positive and negative values. The 
numerical simulations were performed in these two sectors. The analytical considerations 
presented in this section are meant to provide a physical understanding of the main features 
of the early phases of LQC dynamics. Since these do not depend on the sign of 0b, for 
concreteness, in the detailed discussion of analytical considerations we will assume (j)^ to be 
non-negative (which corresponds to numerical simulations reported in Table I) and comment 
on the 0B < case only at the end. 



1. Super-inflation 

Immediately after the bounce, there is a super-infiation phase. Since the potential energy 
at the bounce is very low, this phase of evolution can be well approximated by adding only 
small corrections to the analytically well-understood massless case [12, 15]. From the bounce 
to the end of super-inflation, the pair {H, H) goes from [H — 0, H — 10.28(1 — /^)) to 
{H — 0.93, ij = 0) in Planck units. Thus, this phase is highly dynamical for the Hubble 
parameter. Let us ignore terms of the order Then, the amount of time. At, that 
super-inflation lasts is well approximated by 

A , -^max -^max -^^max /-. . i;2\ n 10 t a 1 r>\ 

At ?a -. fa — = — - fa (1 + f ) fa 0.18 Spi (4.12) 



iiTavg 27r0| 47rpcrit(l - /') 47rp, 



trit 
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How much does the volume change? We can estimate the number of e-foldings from the 
bounce to the end of super-inflation as follows: 

log(iV) = jHdt ^ H,,,At ^ (1 + f) (4.13) 

OTrpcrit 

This implies that the number of e-foldings during super-inflation is only about 1.09. During 
this time the change in the value of the inflaton is given by: 

A</. = /0dt ^ 4vgAt ^ ^ 0.14(l-^)mpi (4.14) 

which is also very small. Thus, although the change in the value of the Hubble parameter 
during super-inflation is dramatic, because the duration of this phase is so short, the total 
changes in the values of the scale factor and the inflaton are also very small. While we 
approximated H and by their average values in these calculations, the final results for 
At and A0 are in excellent agreement with the exact numerical calculations summarized in 
Table I. 

Remark: Dynamics during inflation is often described using an analogy with a damped 
simple harmonic oscillator. This is because when the slow roll conditions hold H is ap- 
proximately constant and equation (4.5) resembles that satisfied by a damped harmonic 
oscillator, with damping parameter ( = 3H/2m. Although the universe does undergo an 
accelerated expansion during super-inflation, since H changes radically, intuition derived 
from a damped harmonic oscillator is not useful in this phase. At the end of super-inflation, 
the Hubble parameter takes its maximum value, i^max = 0.93 mpi. Therefore, the 'friction 
term' C = 3H/2m is extremely large, ( ~ 1.15 x 10^. However, because the kinetic energy 
is very large — approximately half that at the bounce — it still takes an appreciable time for 
the inflaton to slow down. During this rather long period, dynamics is still not mimicked by 
a damped harmonic oscillator. It is only when the potential energy dominates — as is the 
case during the desired slow roll — that the analogy becomes useful. 



2. From the end of super-inflation to turn around 

Let us now consider the post-super-inflation phase. An analytical calculation can be 
carried out in two steps with approximations tailored to each case: i) the phase which com- 
mences at the end of super-inflation and ends when the kinetic energy equals the potential 
energy; and ii) the phase between the time this equality is reached and turn-around of the 
inflaton, i.e., when the kinetic energy reduces to zero. Dynamics has qualitatively difl^erent 
features in these two phases because the flrst is dominated by kinetic energy and the second 
by potential. 

In phase i), the ratio of the potential to the kinetic energy is function a{t) which takes 
values in (0, 1). The approximation consists of replacing it with its 'average' value a^. At 
the end of the calculation, ao is determined by comparing the analytic expression with the 
exact numerical answer in one of the cases given in Table I and then used in other cases. 
We start with the equation b — — 47r70^ from (2.9) and, using the constraint equation (2.8) 
and the approximation, obtain a differential equation containing only b. The solution b{t) 
is given by 

^,/ X N 2.4439 A (t - te) 
cot Xb{t) - cot Xb{te) = ^ (4.15) 

1 ~\~ (Xq 
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TABLE I: Dynamical evolution of (f) and H for various values of / in the case > 0. Analytical 

considerations imply that certain entries should be identically zero; their values are are 0{e.raachine) 
in simulations but denoted by a boldface zero (0) in the Table. Events considered are: the bounce 
point, end of super-inflation (End of SI), equality of potential and kinetic energy {KE = PE), 
turn around = 0), and reaching e = 8 x 10"^. The desired slow roll compatible with the WMAP 
data is not realized if / < 1.17 x 10"^ but is realized for / > 1.25 x 10"^. 

where te is the proper time at the end of super-inflation. It can be substituted back in the 
original Friedmann equation to obtain t — as a function of (f){t) — (p{te)- Finally, we set 
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Method 


= 9.44 X 10"^ 


= 1.02 


= 1.08 


0e = 1.20 


(l)e = 1.33 


0e = 1.46 


Exact Numerics 


4.04 X 10^ 


3.94 X 10^ 


3.87 X 10^ 


3.73 X lO'' 


3.59 X lO'^ 


3.46 X lO'^ 


Analytic Eq(4.17) 


4.00 X 10'' 


3.91 X 10^ 


3.84 X lO'' 


3.71 X 10^ 


3.58 X 10^ 


3.46 X 10^ 



TABLE II: The time interval tsq — teSSSi function of the value (^g of in the infiaton at the end of 
super-inflation: Comparison between the exact numerical result and the analytical approximation. 
As 4>e increases, the difference between kinetic and potential energies at the end of super-inflation 
decreases, whence it less takes less time to achieve equality between potential and kinetic energy. 
The analytical approximation replaced the ratio a{t) between the two energies by an 'average' ag 
during this phase i) after the end of super-inflation. 



t — tEq, the time at which kinetic and potential energies are equal, to obtain: 

5 56 

</'(iEq)-</'(^e) = 0.1628 VrT^(z/Eq-Z/e) wheie sinh(y) = 1 + (— )t. (4.16) 

1 + CKo 

Next, the constraint equation (2.8) determines 0(tEq) in terms of 6(tEq) and we can invert 
the expression (4.15) of b(t) to eliminate it in favor of t^q — te- The result is an expression 
relating tEq — and 0(te) alone. It determines how long the first post- super-inflation phase 
lasts for any given value 4>e of the infiaton at the end of super-inflation: 

(f)e = — , - 0.1628 VI + "0 {VEq-ye) 

» ^-^"^-^"'(^^-o) _ o.l628VTT7^;(s,nh-'(l + 5#(^a>) - 0.88) (4.17) 

^Eq 1 + OCo 

where for simphcity we set te — in the last step. Comparison of this expression with the 
entry / = 1.77 x 10~^ in table Table I, one finds that the fit with the exact numerical result 

is best when ao = 1/3. We will use this value and compare the approximate analytical 
expression with the exact numerical result in all cases covered in Table I. The resulting 
Table II shows that the agreement is very good; the analytical approximation we made by 
replacing a with its 'average' ao works quite well. Table II also shows that the duration 
^Eq — ~ 10^ spi of phasc i) is much longer than that of super-inflation. 

Eq (4.17) determines the duration tEq — te of phase i) as a function of the value of the 
inflation at the end of super-inflation. The same approximation as was used there provides 
expressions for values of H, (p, at the end of phase i) as a function of tEq — t^- 

. l + a„ ., , 0.1152(1 + a,) , 9.5207 x 10^(1 + ^0) 

O.UUl^lEq ~ JeJ ^Eq ~ ^Eq ~ 

(4.18) 

Again, there is good agreement with the exact numerical results. For the Hubble parameter 
a comparison is given in Table III. 



Thus the analytic approximation yields a good portrait of dynamics of all physical quan- 
tities of interest. Note that the Hubble parameter changes very significantly during this 

5q) 



phase also: At the end of super-inflation, we have H ^ 0.93 mpi while at tEq, ~ 10 ^mp\. 



Similarly, while fa 0.9m|i at the end of super-inflation, we have 0(tEq) ~ 4 x 10 ^ 
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Method 


= 9.44 X 10"^ 


0e = 1.02 


= 1.08 


0e = 1.20 


(t>e = 1.33 


</.e = 1.46 


Exact Numerics 


1.02 X 10-5 


1.04 X 10-5 


1.06 X 10-5 


1.10 X 10-5 


1.15 X 10-5 


1.19 X 10-5 


Analytic Eq(4.18) 


1.11 X 10-5 


1.13 X 10-5 


1.15 X 10-5 


1.19 X 10-5 


1.24 X 10-5 


1.28 X 10-5 



TABLE III: Value of the Hubble parameter H at t = tEq as a function of the value Comparison 

between the exact numerical result and the analytical approximation. The analytical approximation 
replaced the ratio a{t) between the two energies by an 'average' Oq during this phase i) after the end 
of super-inflation. The resulting error is less than 10%. Since (pe ^ (pB = /^Amax, H{tEq) increases 
monotonically with /. Since if(tEq) ~ I0~^mp\, general relativity is an excellent approximation 
to LQC at t = tEq- 



Method 


(pe = 9.44 X 10^1 


(pe = 1.02 


(pe = 1.08 


(f>e = 1.20 


(/)e = 1.33 


(/.e = 1.46 


Exact Numerics 


1.26 X 105 


1.24 X 105 


1.21 X 105 


1.19 X 105 


1.15 X 105 


1.12 X 105 


Analytic Eq(4.18) 


1.26 X 105 


1.23 X 105 


1.21 X 105 


1.17 X 105 


1.12 X 105 


1.08 X 105 



TABLE IV: Duration of phase ii), ^ta ~ ^Eq) as a function of (pg: Comparison between the exact 
numerical result and the analytical approximation. The analytical approximation assumes that 
0Avg is given by /? x (p(tEq) where /3 is set equal to 1/4 by agreement with numerical values for 
= 0.944mpi. the analytical expression then reproduces the exact numerical results for other 
values of (pe to within 4%. 



Let us now turn to the subsequent phase ii) which lasts until the inflaton turns around, i.e., 
its kinetic energy vanishes identically. Since the energy density has dramatically decreased 
during phase i) — at tsq, the matter density p(tEq) is of the order IO^^Vpi cases 
given in Table III — general relativity is an excellent approximation in this phase. Also, 
since the kinetic energy is sub-dominant, in contrast to phase i), phase ii) is dynamically 
rather quiescent. At the start of this phase, values of H,(p,(p are given by (4.18). Using 
them and (4.5) we can calculate (p{tEq) and 0(tEq). In all cases summarized in Table III, 
0(^Eq) ~ IQ-^^mli and 0(tEq) ^ 10-^*^mpj. Because of this quiescent behavior we can hope 
to obtain the time duration (txA — ^Eq) until the turn around by a simple calculation: 



TA — tEq J 



A0 _ 4>{tEq)_ 
q) 

^ (4.19) 



0avg /30(^Eq) 
0(^Eq) 



/3{3H{tEq) + m)<P{tEq) P{^H{tEq) + m) 

where ^ is a coefficient (of 0(1)) relating the value oi (p ai t — tEq to its average value in 
phase ii). (As explained in the caption of Table IV, /3 was set to 1/4 in numerical evaluations 
of (4.19).) 

Table IV provides a comparison with the exact numerical values and those calculated 
using (4.19). Note that the quiescent phase ii) lasts even longer than the more dynamic 
phase i). In the case of super-inflation we could use the 'averaging approximation' because, 
although the Hubble parameter is very dynamic, the phase is very short lived. In phase ii) 
we could use it because although the phase is long lived, it is quiescent. Finally, although 
phase ii) lasts ~ 105spi, because H, (p are < 10~^°, < 10"^ t — tEq, the change in H, cp 
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during phase ii) is quite small, again highlighting the quiescent nature of the phase. 



1 




-4 - 



-5 - 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

FIG. 1: Illustrative evolution of the scalar field. In this plot we use a negative value of (pB to 
complement the analytical discussion: (pB = — 5.5mpi and 0b > 0. (See the case / = — 7.35x 10~^ 
in Table V.) The inflaton rolls down the potential quickly at first because the kinetic energy at 
the bounce is large. The kinetic energy equals potential energy at t = tgq = 3.51 x 10'^ spi and 
(j) changes more slowly during the subsequent quiescent phase. The WMAP value e = 8 x 10~'^ 
is reached t = 7.08 x 10^ spi and the end of slow roll, e = 0.1, used in this paper is reached at 
t = 1.24 X lO^spi- Because is negative, in this case the inflaton turns around after the end of 
the desired slow roll, at t = 1.84 x 10^ spy. 



3. Post turn- around 

At the end of phase ii), the inflaton turns around and starts rolling down the potential. 
Since by definition (^(^ta) = 0, the slow roll parameter e also vanishes there. However, 
because if oc 0^, it also vanishes there, whence the slow roll parameter r] need not be small 
and we are not assured the beginning of a slow roll. The desired slow roll, if it is realized in 
the given solution, must occur subsequently, as the inflaton rolls down the potential. 

Recall that, for the desired slow roll to be realized, we must have (f){t{ki,)) ^ 3.15 mpi. This 
point can be reached if and only if ^(tTA) > 3.15 mpi. When is this condition realized and 
when does it fail to be realized in the extreme kinetic energy dominated case under consider- 
ation? The analytical formulas of the last two sub-sections show that 0(tTA) increases with 
/ and the illustrative numerical values we have presented ensure that 0(tTA) > 3.15 mpi 
if 0B ^ 0.937 mpi. (This is a sufficient condition rather than a sharp bound.) Since 
101 G (0, 7.47 mpi) in the extreme kinetic energy dominated case, our result suggests that it 
is likely that the desired conditions are met. However, in the analytical considerations we 
focused only on the 0b > case and moreover made some approximations. Therefore, we 
also carried out a large number of numerical simulations of exact equations to first confirm 
and then sharpen this implication. (Tables I and V contain only a few illustrative results.) 
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FIG. 2: Illustrative evolution of the scale factor for the same case as in Fig 1, i.e., for (j)^ = —5.5 mpi 
and 0B > 0. (See the case / = -7.35 x 10"^ in Table V.) The WMAP value e = 8 x 10^^ reached 
t = 7.08 X 10^ spi and the end of slow roll, e = 0.1, used in this paper is reached at t = 1.24 x 10^ spi. 
The number of e-foldings flattens out soon after this time. 

The simulations were performed in MATLAB using a Runge-Kutta (4,5) ordinary differ- 
ential equation solver (ODE45) with both relative and absolute tolerances set at 3x 10^^^. To 
facilitate computation, the logarithm of the ratio of the volume at a given time to the volume 
at the bounce was chosen to be a fundamental variable. On each solution the preservation 
of the Hamiltonian constraint was verified at the level of the tolerances. For robustness 
sample solutions were also evolved independently with various different tolerances, using 
a set of equations which did not involve the initial volume, and seen to agree with those 
given by the Runge-Kutta algorithm. Time scales for the evolutions were chosen such that 
the numerical noise introduced by the deviation of the Hamiltonian constraint from zero 
would not contribute significantly (i.e. greater than one part in 10^) to the fundamental 
variables. Separate searches were used around the end points of the allowed space, i.e., very 
low values of / for which the dynamical trajectories fail to be consistent with the WMAP 
data. This ensured that the extremely high probability we found for the LQC dynamics to 
be compatible with the WMAP data is in fact a lower bound. 

As discussed in section IV A, for the quadratic potential the WMAP data implies that 
there was a time t{k^,) in the early history of the universe when e = 8 x 10~^ and = 3.15 mpi 
within error bars (of < 4.5%). In LQC we can specify the initial data at the bounce. 
Numerical simulations let us answer the following question for the extreme kinetic energy 
dominated bounce: What are restrictions on the initial data at the bounce for the ensuing 
LQC dynamical trajectory to meet the WMAP constraint? Full numerical simulations bear 
out conclusions suggested by Tables 1 and V: 



Thus, in the extreme kinetic dominated case, we are guaranteed that the event E — the 
desired slow roll compatible with WMAP observations — will be realized in the solution under 



WMAP constraint is met provided 




(4.20) 
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/ 


Event 


<^ 




H 


H 


t 


e 


—7.02 X 10 ° 


Bounce 
End of SI 
KE = PE 

= 
e = 8 X 10-3 


-5.25 
-5.11 
-3.15 
5.42 X 10-2 
NA 


9.05 X 10-^ 
6.39 X 10-1 
3.80 X 10"^ 

NA 




9.26 X 10-1 
1.10 X 10-^ 
2.25 X 10-^ 
NA 


1.03 X lOi 


-1.82 X 10-1° 

NA 




1.80 X 10-1 

3.81 X 10^ 
1.71 X 10^ 

NA 


oo 


1.50 


NA 


—7.22 X 10 ° 


Bounce 
End of SI 

KE = PE 

4> = o 

e = 8 X 10-3 


—5.40 
-5.26 

-3.31 
5.64 X 10~2 
NA 


9.05 X 10-1 
6.39 X 10-1 
4.00 X 10-^ 


NA 




9.26 X 10-1 
1.16 X 10-5 
2.14 X lO-'^ 
NA 


1.03 X lOi 


-2.01 X 10-1° 


NA 




1.80 X 10-1 
3.62 X 10^ 
1.79 X lO'^ 
NA 


oo 


1.50 


NA 


-7.35 X 10 


Bounce 
End of SI 
KE = PE 

(j) = 
e = 8 X 10-3 


-5.50 
-5.36 
-3.41 
5.51 X 10-2 
-3.14 


9.05 X 10 ^ 
6.39 X 10-1 
4.12 X 10-6 


1.97 X 10-'^ 




9.26 X 10-1 
1.19 X 10-5 
2.21 X lO-'^ 
7.80 X 10-6 


1.03 X lOi 


-2.14 X 10-1° 


-4.86 X 10-13 




1.80 X 10-1 
3.51 X 10^ 
1.84 X lO'^ 
7.08 X 105 


oo 

1.50 


8 X 10-3 


—7.69 X 10 


Bounce 
End of SI 
KE = PE 

(/) = 
e = 8 X 10-3 


-5.75 
-5.61 
-3.68 
5.71 X 10^2 
-3.14 


9.05 X IQ-i 
6.39 X 10-1 
4.44 X 10-6 


1.97 X 10-'^ 




9.26 X IQ-i 
1.29 X 10-5 
2.11 X lO-'^ 
7.80 X 10-6 


1.03 X lOi 


-2.48 X 10-1° 


-4.87 X 10-13 




1.80 X 10-1 
3.26 X lO"! 
1.97 X lO'^ 
2.03 X 106 


oo 

1.50 


8 X 10-3 


—7.92 X 10 ° 


Bounce 
End of SI 

KE = PE 

= 
e = 8 X 10^3 


—5.92 
-5.78 

-3.85 
5.48 X 10-2 
-3.14 


9.05 X 10-1 
6.39 X 10-1 
4.66 X 10-6 


1.97 X lO-'^ 




9.26 X 10-1 
1.35 X 10-5 
2.22 X 10-^ 
7.80 X 10-6 


1.03 X lOi 


-2.73 X 10-1° 


-4.86 X 10-13 




1.80 X 10-1 

3.10 X 10^ 
2.07 X lO'^ 
2.94 X 106 


oo 

1.50 


8 X 10-3 


—8.02 X 10 ° 


Bounce 
End of SI 
KE = PE 

= 
e = 8 X 10-3 


-6.00 
-5.86 
-3.94 
5.51 X 10-2 
-3.14 


9.05 X 10-1 
6.39 X 10-1 
4.75 X 10-6 


1.97 X 10-'^ 




9.26 X 10-1 
1.38 X 10-5 
2.21 X 10-'^ 
7.80 X 10-6 


1.03 X lOi 


-2.84 X 10-1° 


-4.86 X 10-13 




1.80 X 10-1 
3.04 X 10^ 
2.11 X 10'^ 
3.36 X 106 


oo 

1.50 


8 X 10-3 







TABLE V: Dynamical evolution of and H for various values of / in the case 0b < 0. Analytical 

considerations imply that certain entries should be identically zero; their values arc arc 0{emachine) 
in simulations but denoted by a boldface zero (0) in the Table. Events considered are: the bounce 
point, end of super-inflation (End of SI), equality of potential and kinetic energy {KE = PE), 
turn around (0 = 0), and reaching e = 8 x 10-3. Because 0b < 0, the turn around happens after 
the desired inflation ends. The desired slow roll compatible with the WMAP data is not realized 
if I/I < 7.22 X 10-6 but is realized if |/| > 7.35 x 10-6. 
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consideration if f is outside the interval (—7.35 x 10^®, 1.25 x 10~^). (Again, this is only 
a sufficient condition.) Since full range of / in the kinetic energy dominated case under 
consideration is (—10^^, 10^^), this result suggests that E will be realized in 'more than 
half the trajectories' in this case. We will sharpen this statement using the normalized 
Liouvillc measure in section IV D. 

The precise numbers in the sufficient condition (4.20) depend on the sign of (pB at the 
bounce. We will conclude with a brief discussion of the origin of the difference between 
the > and 0b < cases. The origin of the asymmetry hes in the fact that (because 
of symmetries of field equations) we have restricted ourselves to data with 0b > 0. So, if 
0B < 0, the inflaton is already rolling down the potential at the bounce. If initially it is 
sufficiently high up in the potential with, |0b| > 5.5 mpi, friction can slow it down sufficiently 
for the desired slow roll to begin at ~ — 3.15mpi. But if \4>b\ < 5.5mpi, then the kinetic 
energy is too large for the friction term to slow it sufficiently for the desired slow roll to 
commence before it reaches the bottom of the potential. Then it starts climbing up on the 
(f) > branch but does not acquire a value higher that ~ 3.15 required for the dynamical 
trajectory to pass through the region satisfying the WMAP constraints. In the 0b > case, 
by contrast, if 0b ^ 0.94, the kinetic energy at the bounce is sufficient to propel the inflaton 
to values higher than 3.15mpi on the 0b > branch of the potential to allow for the desired 
slow roll. 

In subsection IV C 2, we will show that the likelihood of attaining the desired slow roll is 
in fact 1 if the bounce is not extreme kinetic energy dominated and in section IV D we will 
calculate the precise probability for the occurrence of E on the entire space S of solutions 
using the normalized Liouville measure. 

C. Generic LQC solutions and the WMAP slow roll 

In section IV CI we continue our description of qualitative features of LQC dynamics, 
now for bounces with / > 10~^. This discussion is rather sketchy compared to that of 
section IV B because in this case there are no subtleties with respect to the WMAP data: as 
we show in section IV C 2, all dynamical trajectories meeting this condition at the bounce 
satisfy the WMAP constraints. 

1. Evolution of other initial data: Qualitative features 

We will highlight only those features which are qualitatively different from the kinetic 
energy dominated bounce discussed in detail in section IV B. Again, for concreteness, we 
will focus on the case when 0b and 0b are both positive. Let us begin with the intermediate 
case 10~^ < / ^ 1/a/2 where the kinetic energy still exceeds the potential energy at the 
bounce but does not dominate it as in section IV B. The super- inflation era is similar to that 
described in section IV B 1 but because 0b is now lower, the phase lasts longer. The inflaton 
climbs up the potential but change in its value is again small. The Hubble parameter, on 
contrast, is again dynamical. However, the post super-inflation dynamics exhibits signiflcant 
differences. For, now the value of 0b is higher and 0b lower while, as before, H assumes its 
largest value at the end of super- inflation. Therefore, the coefficient of friction, ( — 3H/2m, 
is again large but there is less kinetic energy to lose before reaching the turn-around., which 
is now reached within 10— lOOspi after the bounce. Consequently, now the change 0(iTA) —0b 
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is negligible, a key feature not shared by regime (i). 

Next, let us examine the / > 1/ \/2 where the potential energy is greater than the kinetic 
energy. Now, the LQC effects dominate for an even longer time. Again, because > 0, the 
inflaton climbs up the potential but turns around earlier and earlier and the super-inflation 
phase lasts longer and longer as / increases. Now the turn around {(f) = 0) will occur during 
super-inflation! The change (0(tTA) ~4>b) is even more negligible because the kinetic energy 
at the bounce is lower than that in the / < 0.5 case. The slow roll conditions are easily 
met soon after turn-around. A difference from this phase of slow roll and that for / > 1/ \/2 
is that H continues to grow during the slow roll because we are still in the super- inflation 
phase. There is an enormous number of slow roll e-foldings already in the deep Planck 
regime where the matter density is greater than half the critical density. 

2. Compatibility with the WMAP data 

In the extreme kinetic energy domination we could carry out numerical simulations start- 
ing from the bounce until the end of the desired slow roll. For / > 0.1, this becomes quite 
difficult even with the truncation error as small as one part in 10^^ because the super-inflation 
phase can become so long that the required evolution spans ~ lO^^spi or more. Therefore, 
to maintain the precision we had in the kinetic dominated case, another line of reasoning is 
necessary. 

X 10"' 



8.28 - 
8.26 - 




8.22 - 
8.2 - 
8.18- 

I I \ I I \ I I 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

df dt I (^=-5.5 

FIG. 3: WMAP data implies that, for a quadratic potential, at the time t{k-^) at which the reference 
mode ki, exits the Hubble horizon, \<p\ = 3.15 mpl and e = 8 x 10~^ within error bars of ~ 4.5%. We 
wish to analyze whether generic dynamical trajectories pass through this small neighborhood in 
the phase space. For reasons discussed in section IV C 2, we consider trajectories with initial data 
with a fixed, low value of (/) (which we take to be (/> = — 5.5mpi), allowing all permissible values of 
<j). The plot shows that in all the resulting dynamical trajectories, at the time when (p = — 3.15mpi, 
we have e ~ 8 x 10^^ with fluctuations less than 2 x 10~^. This implies that all LQC trajectories 
with I/I > 7.35 X 10~^ pass through the desired, small region of the phase space. 

Motivated by our result (4.20), in this subsection we consider initial data at the bounce 
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with |0b| ^ 5.5mpi7 We will first argue that the resulting dynamical trajectory must pass 
through the point with |0| = 5.5 mpi with some kinetic energy. Since this is a low value of 
1 01, the inflation will reach the WMAP value = 3.15mpi rather quickly. So it is feasible 
to carry out a numerical simulation for dynamics between |0| = 5.5 mpi and \(f>\ — 3.15 mpi 
with the desired, high degree of accuracy. The question for numerics is then: Is the value 
of e at the time when |0| = 3.15 mpi close to the WMAP value e = 8 x 10~^ within the 
4.5% error bars? For initial data at the bounce for which |0b| ^ 5.5mpi, the super-inflation 
phase will last very long time and, had we begun the simulation at the bounce, we would 
have lost accuracy by the time the inflaton reaches |0| = 3.15 mpi. We avoid this problem 
by giving an analytic argument that the inflation must reach |(/)| = 5.5 mpi and carry out 
numerical calculations for dynamics only between = 5.5 mpi and |0| = 3.5 mpi. 

Let us begin with the negative branch of the 0b space, i.e., flrst consider the case when 
— 5.5mpi > (f)B> — |0max| = —7.47 X lO^mpi. In this case, immediately after the bounce the 
inflaton will roll down the potential. Intuitively, it may seem obvious that in this descent it 
must encounter the point (f) = — 5.5mpi with some kinetic energy. But a priori there are two 
possibilities that may prevent this occurrence. First, we have to examine the possibility that 
the inflaton may come to rest for some value (j) = (po < —5.5 mpi and just stops there. But 
the equation (4.5) satisfled by the inflaton imphes that at the instant = 0, (f) — — m^0o > 
(since 0o 7^ 0). Therefore the inflaton cannot just stay at 0o < 0. A more subtle possibility 
is that the inflaton asymptotically approaches 0o < — 5.5mpi, with approaching zero, but 
never actually reaches it in a flnite time. That is, the limiting value of could be zero at 
= 00 but this could happen at time to — 00. But by integrating (4.5) between the bounce 
time tB and the hypothetical time to at which is to vanish, it is easy to show that: 

to-tB< (4.21) 

m^|0B| 

so that to is necessarily finite. Thus, if 0b < — 5.5mpi (and by assumption 0b > 0), the 
infiaton must slide down the potential and reach the value — — 5.5mpi with some kinetic 
energy. 

Next, consider the case when 0b > 0. In this case, our focus will be on initial data at 
the bounce with 0b > 5.5 mpi. Since we again have 0b > 0, initially the inflaton now climbs 
up the potential. However, it again follows from (4.5) that it reaches = in a flnite time 
to and at that time = —m'^(f) < (since is now positive). So this is the turn around 
point and the inflaton rolls down the potential. Again, by the same reasoning as before, as 
it rolls down, the inflaton cannot come to rest for a value 7^ 0. Therefore as it rolls down, 
it must pass through the point = 5.5 mpi, for some value of 0. To summarize, dynamical 
trajectories for every initial data at the bounce with |0b| > 5.5 mpi eventually encounter a 
point at which |0b| = 5.5 mpi. 

We used high precision numerics to analyze the LQC dynamics following this event. The 
value of at this event can be arbitrary. We sampled the full range, (0 = 0, = 0.90 m|j), 
first using uniformly distributed thousand data points and then logarithmically distributed 
thousand points. Using = — 5.5mpi and each of these values of as initial data solved the 
full set of LQC equations numerically. The key questions then are: i) Do all these dynamical 



This space has an overlap with the / < 10 ^ case already considered in detail in section IV B. This 
ensures that there are no 'gaps', i.e., we will cover the full space of initial data at the bounce. 
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trajectories eventually pass through a phase space point at which = — 3.15mpi; and, ii) 
Within the WMAP error bars, which of them have e = 8 x 10~^ at the time when assumes 
the value 3.15 mpi? As one would expect from the analytical considerations discussed above, 
the answer to the first question is in the affirmative. As the infiaton rolls down the potential 
from = — 5.5mpi, it necessarily encounters the value = — S.lSmpi. The answer to the 
second question is summarized in Fig. 3: in each of these two sets of 1000 simulations, 
at the time the infiaton assumes the value = — 3.15mpi, the slow roll parameter e takes 
values in the range 

e = 8 X 10"^ ± 2 X 10~^ . (4.22) 

That is, each of these trajectories passes through the small portion of the phase space 
compatible with the WMAP data. 

Finally, let us consider the complementary case with initial data at the bounce satisfying 
(pB > 5.5 mpi. In this discussed above, the infiaton first rises up the potential and 

then rolls down. Therefore now the numerical simulation should start with = 5.5 mpi and 
4> e (— O.QOmpi, 0). But because of the symmetry on the space of solutions noted in the 
beginning of section IV B, these solutions can be obtained just by reversing the sign of the 
solutions (j){t) obtained from numerical evolutions starting from the initial data (f) = —5.5 mpi 
and 4> G (0, 0.90 mpi). Thus, all the dynamical trajectories are compatible with the WMAP 
data also in this case. 

To summarize, combining the results of sections IV B 3 and numerical simulations just 

discussed, we can conclude that the LQC dynamical trajectories resulting from any initial 
data at the bounce surface with |0b| > 5.5 mpi realizes the constraints on values of fields at 
t = t{ki,) imposed by the WMAP data. Thus, it is only when the bounce is in the extreme 
kinetic dominated regime that some of the dynamical trajectory — namely those that violate 
(4.20) can fail to meet the WMAP constraint. 



D. Probability of the desired slow roll in LQC 

We are now ready to combine results of sections IIIB, IV B and IV C to calculate the a 
priori probability of realizing the desired slow roll in LQC. 

In the terminology of section III B, the event E of interest is the passage of the dynamical 
trajectory through the small region in the phase space singled out by the WMAP data and 
our task is to find the relative volume of the region R{E) of the space of solutions in which 
E occurs. The calculation requires a normalized measure. In section IIIB we found that, 
thanks to a suitable gauge fixing, each equivalence class of physically distinct LQC solutions 
can be characterized completely by the value 0b of the infiaton at the bounce surface. Since 
h\ = n/2 at the bounce, (3.19) implies that the total measure on the space §o of physically 
distinct solutions is given by 

r /•"/'max 0„ 1 0_ 

where we have used the fact that we have a quadratic potential and where, for definiteness 
we have set the volume Vo at the bounce to 1. (Recall that Vg simply rescales the measure 
by a constant and therefore drops out in the calculation of probabilities.) The probability 
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P{E) of the occurrence of our event E is then given by (3.20): 

^(^) = ^/ (^-4AW)M0 (4.24) 

-'^ Jx{E) ^ 

where X[E) is the interval on the 0B-axis corresponding to the physically distinct LQC 
solutions in which E occurs. (For details, see section IllB.) 

In sections IV B and IV C we found that a sufficient condition for the event E to occur 
is that the initial data of the solution at the bounce should satisfy 

0B ^ [-5.5 mpi, .94mpi] or, cquivalently, / ^ [-7.35 x 10"^ 1.25x10"^] (4.25) 

Therefore the probability that E does not occur is bounded by the length of the interval 
[— 5.5mpi, .94mpi] in the 0B-axis w.r.t. the measure given by the 1-form a): 

1 r4'+ r4'- o ^ 

P(E not realized) < ^(/ +/ ) - 47rV^V^) ' # 

^ Jo Jo ^ 

< -{Ujl-fl + sm-'f+) + -ifJl-fl + sm-'f.) 

< 2.74 X 10"^ (4.26) 

where 0+ = .94mpi, 0_ = 5.5 mpi and f± = 0±/0max- Thus, the probability that the desired 
slow roll does not occur in an LQC solution is less than three parts in a million. Hence a 
great deal of fine tuning would be necessary to avoid the slow roll inflation that meets the 
WMAP constraints. 

For simplicity, throughout our analysis, we set the cosmological constant A to be zero. 
But it is completely straightforward to include it since it just shifts the zero of our potential. 
Since the sign of the observed cosmological constant is positive, for any given value of 0, its 
inclusion would have the effect of increasing the fraction of the total energy density in the 
potential. This in turn would lower the vahics of |(/)b| that are needed in (4.20) to ensure 
that the trajectory meets the WMAP constraint. Thus the probability of achieving the 
desired slow roll would further increase. However, since the observed value of A is so small, 
this effect will not be noticeable even at the level of numerical accuracy used in this paper. 
We also focused on the spatially flat FLRW models because they are observationally favored. 
But it is not difficult to include spatial curvature. Again, given observational constraints, 
we expect that this inclusion will not significantly alter our main conclusions [34]. In these 
respects, the results are robust. 

What about the choice of the inflaton potential? All our detailed considerations hold only 
for a quadratic potential. Recall however that a generic potential V{(f)) is well-approximated 
by a quadratic one near its minimum. Suppose V{(j)) has a single minimum (po and that 
|V(0) - V{(j)o) - (l/2)mV^| -C V{(j)) - vl^o) for aU e [-5.5 mpi, 5.5 mpi] where m ^ 
1.21 X 10~^mpi. Then our considerations will apply: As the inflaton shdes down from = 
5.5mpi it will enter a phase of slow roll inflation compatible with the WMAP constraints. 
As an illustration, let us suppose that the potential is a combination of a quadratic and a 
quartic parts: 

F(0) = imV' + (4.27) 

Then, our assumption is met if A <^ 9.68 x 10~^"^. For a rough comparison, note that the 
WMAP bound for a pure quartic potential is 4.4 x 10~^^. Thus, it is not unreasonable to 
expect that, at a qualitative level, our conclusions on probabilities will continue to hold for 
a wide class of physically interesting potentials. 
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V. DISCUSSION 

In the main body of this paper we analyzed inflation in the observationally most interest- 
ing case of the k=0 FLRW cosmologies using the framework of LQC. Our goal was two-fold. 
First, assuming a quadratic potential, we examined in detail the effective LQC dynamics in 
presence of an inflaton with standard potentials. Second, we used this dynamics to calculate 
the a priori probability of realizing the desired slow roll inflation that is compatible with 
the WMAP data. 

Let us begin with dynamics. Since the big bang singularity is replaced by the big bounce, 
and since all physical fields are regular at the bounce, we could specify initial conditions 
at the bounce surface and explore the resulting dynamics. We allowed all possible initial 
data, subject only to the Hamiltonian constraint. We found that the subsequent dynamics 
is sensitive to the distribution of energy density at the bounce between kinetic and potential 
parts. In all cases, there is a universal super-inflation phase that follows the bounce; indeed 
it persists even when there is no potential at all. For concreteness, let us focus on the case 
when 4> and 4> are both positive at the bounce so the inflaton is climbing up the potential. 

We first discussed in detail the case in which less than one part in 10^° of the total 
energy density at the bounce is in the potential because this is the only case in which some 
of the solutions can fail to exhibit the desired slow roll. In this case, the super-inflation 
phase is very short lived. However, it is also extremely dynamic: during the fraction of a 
Planck second of super-inflation, the Hubble parameter increases dramatically from zero to 
its maximum allowed value, 0.93 mpi. Because H is so large at the end of super-inflation, the 
friction term in the equation of motion of the inflaton is also large whence the inflaton loses 
kinetic energy as it climbs up the potential. After ~ 10^ spi, the potential energy equals 
kinetic energy and then dominates it (where, as before, spi denotes Planck seconds). After 
another ~ 10^ — 10^ spi, vanishes, the inflaton turns around and starts climbing down 
the potential. In a 'majority' of solutions — but not all — it soon enters the desired slow roll 
that is compatible with the 7- year WMAP observations [6]. 

In the intermediate case, the kinetic energy is still greater than the potential at the bounce 
but does not overwhelm it. Now dynamics undergo the same qualitative phases. However, 
the super-inflation phase lasts a bit longer and the turn around occurs much sooner, only 
about 10 — lOOspi after the bounce. This is because the initial kinetic energy is smaller than 
that in the first case while the strength of the friction term is the same because the Hubble 
parameter again takes its maximum value at the end of super-inflation. Because the turn 
around occurs much sooner, in contrast to the situation in the extreme kinetic domination, 
the increase in the value of (p between the bounce and the turn- around is very small. Finally, 
if the potential energy density exceeds the kinetic energy density at the bounce, dynamics 
is very different. Now the LQC effects dominate in the sense that the super-inflation phase 
lasts longer and its duration increases rapidly as the fraction of total energy in the potential 
increases. Furthermore, the turn around occurs already during super-inflation! Our first 
goal was to present these qualitative differences in the pre-slow-roll dynamics. 

The WMAP data puts very strong constraints on the values of fields at the time t{k^) 
when the mode with wave number k^, exits the Hubble radius during the desired phase of 
inflation. The key question is: Are 'most' dynamical trajectories such that these constraints 
are met at some time during the post-bounce evolution? We showed that the answer is in 
the affirmative unless the bounce is kinetic energy dominated as in section IV B. Even in 
this case, the WMAP constraints are met unless fraction of the total energy at the bounce 
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which is in the potential is less than 7.35 X 10~^. In this precise sense the tiny region of phase 
space selected by the WMAP data serves as an attractor to LQC dynamics. This may not 
seem surprising because inflationary trajectories are known to be attractors also in general 
relativity. However, in LQC there are strong quantum gravity efi^ects in the super-inflation 
phase as well as in the phase that follows immediately after super-inflation. These could well 
have spoiled the attractor behavior and a large fraction of the data specified at the bounce 
could well have given rise to trajectories that steer away from the WMAP region. That this 
does not happen is noteworthy. 

The attractor behavior is a qualitative feature of dynamics. It does not provide a sharp 
quantitative estimate on the likelihood of realizing the desired slow roll. To obtain this 
estimate, one needs a measure on the space of (physically distinct) solutions, such that the 
total volume of this space is finite. As explained in section I, a natural strategy [7-9] is to 
use the Liouville measure on the phase space. This leads to a priori probabilities of events, 
such as the occurrence of the desired slow roll. However, volume of the 2-dimensional space 
S of sohitions with respect to the resulting d/tL is infinite. But this infinity is physically 
spurious: It occurs because S is acted upon by a gauge group Q and the length of 1- 
dimensional orbits of Q is infinite. For a large class of potentials the space §>/Q of physically 
distinct solutions is in fact a bounded interval of the real line. Furthermore there is a natural 
1-form Ua on S (constructed from the symplectic structure and the gauge vector field G") 
which is everywhere orthogonal to G". So it is a natural candidate to induce a volume 
element on the 1-dimensional space S/Q. Unfortunately, it does not project down to S/Q 
(because Cc^a = 0). Therefore one is led, instead, to 'fix a gauge', i.e., hft S/^ to 

a suitable 1-dimensional cross-section in § and carry out integration there. The problem in 
general relativity is that there is no natural family of lifts and the volume element on S/^ 
depends heavily on one's choice. Thus, there is an inherent ambiguity in the calculation 
of probabilities. Thanks to the existence of the bounce surface, in LQC one does have a 
natural family of lifts and the probability of occurrence of physical events is independent of 
the choice.^ We used this choice to calculate the a priori probability of realizing the desired 
slow roll, compatible with the WMAP data. We found that the probability is greaier than 

0. 999997 in LQC. As emphasized in section I, these are just 'bare probabilities' and better 
estimates of this occurrence will require astutely chosen physical inputs. This will require 
not only a more complete theory [3], but a deep understanding of that theory to separate 
essential inputs from non-essential ones. However, since the 'bare probability' is so close to 

1, it is a huge burden on any theory to come up with new inputs that significantly change 
the answer. 

We emphasize, however, that our results should not be interpreted to mean that suffi- 
ciently long slow roll infiation is inevitable in LQC: it is inevitable only under the additional 
assumption that there is a phase in which matter density is dominated by that of a scalar 
field in a suitable potential. So far LQC has not provided a mechanism to create either the 
scalar field or the potential. Although there have been intriguing suggestions that this may 
naturally occur if one promotes the Barbero Immirzi parameter 7 of loop quantum gravity 



In this construction, the analog of the bounce surface in general relativity would be the singularity itself. 
It is far from obvious how to work at the singularity although ideas proposed in [35] may provide a 
natural mathematical construction. However, a priori, the physical meaning of such a construction would 
be obscure because one would have to assumes that Einstein's equations, without quantum corrections, 
are valid all the way to the singularity. 
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to a dynamical field [36], very substantial work is needed to have confidence that they can 
be transformed into a concrete, internally complete and viable scenario. 

The primary purpose of our analysis was to understand what LQC has to say about 
inflationary scenarios in the early universe. But it is instructive to compare and contrast 
the predictions of LQC with those of general relativity, discussed in the literature. This 
leads us to several distinct points. First, there is some controversy as to which 'problems' 
inflation solves and which it does not [1-3]. Our analysis does not shed any new light 
on these issues. Rather, as mentioned above, our focus is on the hkehhood of the desired 
slow roll inflation. Second, while some of the literature [4, 7-9] uses the Liouville measure 
on the space of solutions to compute the a priori probability, in [2] (and in the earlier 
literature cited therein) a different measure is used. The likelihood, of course, can and does 
depend sensitively on the choice of the measure. In addition, since the choice made in [2] 
is not preserved by the Hamiltonian flow, a choice of time slice is made. Therefore, there is 
some discussion in the literature on the dependence of the final results on the choice of the 
measure and the time slice used in the evaluation of probabilities (see, e.g., [1, 3, 4]). The 
Liouville measure is 'canonical' and, furthermore, preserved under time evolution. However, 
as explained above, the total measure on the space of solutions § is infinite and, to extract 
meaningful probabilities, one must introduce an additional structure, a suitable lift of E>/Q 
into S. Conceptually, this is equivalent to fixing a 'time gauge,' e.g. by working with the 
initial data of the solution at the instant at which the matter density is a given constant, 
Po- In this sense, there is a common limitation. In [2] (and in earlier works dating back 
to [32]) one uses the time shce at Planck time (but still works with equations of general 
relativity), while in [4] it is argued that one should use a slice corresponding to a much 
later time. Even if one were to decide to use the Liouville measure, a priori probabilities 
do depend sensitively on this choice because of the nature of infiationary dynamics [10]. 
Prom our perspective, since there is no 'canonical' choice of time slice in general relativity, 
calculations of probabihties have an inherent ambiguity.^ In LQC, by contrast, since the 
bounce surface provides a canonical 'time slice' the ambiguity can be naturally resolved. 
Finally, recently there have been several interesting discussions of infiationary scenarios in 
LQC (see, e.g., [37-42]). However, the primary focus of all but one of these papers is 
on phenomenological and observational issues rather than on measures and calculations of 
probabilities. The one exception is [37]. However, in that reference, super-infiation (as well 
as the bounce) was ignored. 

We will conclude with an observation on dynamics immediately following the bounce. In 
the standard infiationary scenarios, matter density is some 11-12 orders of magnitude smaller 
than ppi at the onset of infiation and from the quantum equations of full LQC we know 
that the use of quantum field theory in curved space-times is fully justified in this regime. 
However, there are several conceptually important issues which require an understanding of 
dynamics before this era is reached. Perhaps the most important among them is the issue 
of the initial state of quantum fields representing perturbations. As explained in section I, 
currently the state is simply postulated to be the Bunch-Davis vacuum for the relevant modes 
(which have co-moving wave numbers in the range {kg, ~ 200/Co))- While this assumption 



There is another difference between references [2] and [4]: while [2] focuses on the (phenomenologically 
favored) k=0 model, [4] focuses on the (spatially closed) k=l model. However, since the strategy used 
in [4] to 'regularize' the Liuoville measure is strongly motivated by the gauge considerations in the k=0 
case, this difference is less significant for the drastic disparity in the final results. 
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is physically motivated, it is nonetheless important to arrive at this state starting from 
more fundamental initial conditions. In general relativity, these conditions would have to be 
specified on the singularity and furthermore we know that we cannot reliably use equations 
of general relativity in the Planck regime. LQC on the other hand is based on loop quantum 
gravity, a candidate theory of full quantum gravity. Therefore, not only is there a basis to 
trust its quantum equations but, as we saw in section II, they have already provided us a 
wealth of new information on physics at the Planck scale. Therefore, there is good motivation 
to use quantum field theory on the cosmological quantum space-times of LQC [43] to study 
the initial conditions and evolution of quantum fields representing perturbations from the big 
bounce until the onset of inflation. Is there perhaps a natural initial condition we can impose 
on the quantum state of these perturbations at the bounce using e.g., ideas developed in [44] 
and/or exploiting the fact that the Hubble parameter vanishes there? What would such an 
initial state evolve to at the onset of slow roll? Would it be sufficiently close to the Bunch- 
Davis vacuum (for the modes of interest) to be phenomenologically viable, or, would it be so 
different that it is already ruled out observationally? If it close, is it perhaps too close to be 
observationally indistinguishable from the Bunch-Davis vacuum or is it sufficiently different 
to lead to an observational test of LQC? These questions are being currently analyzed using 
the detailed dynamics of the LQC inflationary space-times presented in section IV B [17]. 

Finally, because almost all LQC solutions are compatible with the WMAP data, one 
might first think that the chances of constraining quantum gravity from observations of 
the early universe are very small. However this is not correct: to make contact with these 
observations one needs not only the 'background' space-time but also quantum fields repre- 
senting perturbations off this background. Considerations of the previous paragraph suggest 
a concrete direction to confront quantum field theory on the LQC quantum space-times with 
observations. Such a confrontation could well lead to interesting, testable predictions of loop 
quantum gravity and/or constraints on this theory. 
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